Reexamination of helium fine structure 
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In order to explain discrepancies between theoretical predictions and experimental data for the helium fine 
structure, we check and recalculate all theoretical contributions up to orders m a 7 and m 2 /M a 6 . The previous 
result for the m a 7 correction is improved by a much more accurate calculation of relativistic corrections to the 
Bethe logarithm. The theoretical values of the 2 3 Po — 2 3 Pi and 2 3 Pi — 2 3 P2 fine structure intervals in helium 
are, correspondingly, v 01 = 29 616 946.2(1.6) kHz and v 12 = 2 291 177.3(1.6) kHz, with the uncertanties 
being due to higher-order effects. For the small interval v\%, the theoretical value agrees with the experimental 
data, whereas for the large interval vqi, a discrepancy of about 3 standard deviations is present. 
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I. INTRODUCTION 

The fine structure splitting of the 2 3 P level in helium has 
long been an attractive subject of theoretical and experimen- 
tal studies. One of the reasons for this interest is that the fine 
structure, being an intrinsically relativistic effect, is propor- 
tional to a 2 Ry and thus provides an opportunity for the deter- 
mination of the fine structure constant a from a comparison 
of theoretical predictions with experimental data. A series of 
measurements of the helium fine structure has been performed 
during the last decade [1-7], with the current accuracy being 
on the level of 25 ppb. For theory, to reach an adequate level 
of precision in a description of a three-body system is a chal- 
lenging problem. 

Despite considerable calculational efforts of last years to 
provide an accurate theoretical determination of the fine struc- 
ture of helium, the current status of theory can hardly be con- 
sidered as satisfactory. Recent calculations [8, 9] demon- 
strated a significant discrepancy with the experimental data, 
the difference for the large (small) fine structure interval being 
about 10 (6) times larger than the total nonlogarithmic contri- 
bution to order ma 7 . It seems unlikely that such difference 
can be explained only by higher-order effects. 

The theory of the helium fine structure up to order m a e has 
been confirmed by at least two independent calculations and 
thus can be considered as established. The logarithmic part of 
the m a 7 contribution has also been calculated independently. 
The only corrections that are not yet checked by different eval- 
uations are the recoil contribution to order m 2 /M a 6 and the 
nonlogarithmic correction to order ma 7 . These corrections 
will be the main subject of the present investigation. 

The theoretical description of the fine structure to order 
m a 7 is a difficult task. Within the logarithmic accuracy, this 
was first done by Zhang et al. [10] and later confirmed by 
one of the authors (K.P.) [11]. An important part of the non- 
logarithmic correction to order m a 7 was calculated by K.P. 
and Sapirstein [12]. The problem of derivation of the com- 
plete set of effective spin-dependent operators to order m a 7 
was addressed by Zhang in a series of works [13-15] within 
the equal-time variant of the Bethe-Salpeter formalism. The 
derivation based on the dimensionally regularized NRQED 
was reported recently by K.P. [9], who noted several mistakes 



and inconsistencies in the previous derivation by Zhang. 

In the present investigation, we give a detailed account of 
the derivation of the total contribution to order m a 7 , first re- 
ported in Ref. [9], and present a recalculation of all corrections 
up to orders m 2 /M a 6 and m a 7 . Particularly, we perform an 
evaluation of the relativistic correction to the Bethe logarithm, 
which improves significantly upon the first calculation in Ref. 
[12]. 

The paper is organized as follows. In Sec. II and III we give 
a short summary of the general formulas for the helium fine 
structure up to order ma e . The derivation of the ma 7 cor- 
rection is presented in Sec. IV. The next section describes the 
numerical approach and reports the numerical results. Sec. VI 
contains the summary of all contributions to the helium fine 
structure and the discussion of the present status of theory and 
experiment. 

The relativistic units are used in this paper, h = c = eo = 1 
and e 2 = Ana. 



II. LEADING-ORDER FINE STRUCTURE 

The dominant contribution to the helium fine structure is 
induced by the spin dependent part of the Breit-Pauli Hamil- 
tonian, which is, for an infinitely heavy nucleus, 
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where r = f\ — r% and we have included the effects of the 
electron anomalous magnetic moment (amm) a e , 
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Expanding the amm prefactors in Eq. (1), Hf s can be written 
as a sum of operators contributing to different orders in a, 



These operators were later re-derived in a more simple way 
using the effective field theory in Refs. [11, 13]. The result is 
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Here and yield the complete fine-structure con- 
tributions of order ma A and ma 5 , respectively, whereas 

#fe 6 Lm and ^felmm ^ ^ amm P 3 ^ of the correspond- 
ing higher-order operators. 

The leading effect of the finite nuclear mass is conveniently 
divided into three parts, termed as the mass scaling, the mass 
polarization, and the recoil operators. The effect of the mass 
scaling is accounted for by including the prefactor (m r /m) 3 
into the operator _fff s , where m r is the reduced mass for the 
electron-nucleus system. The effect of the mass polarization 
can be accounted for to all orders by evaluating expectation 
values of all operators on the eigenfunctions of the Shrodinger 
Hamiltonian with the mass-polarization operator (m r /M) pi ■ 
P2 included. The third effect is induced by the recoil addition 
to the Breit-Pauli Hamiltonian, 
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IH. ma 6 CONTRIBUTION 

The m a e contribution to the helium fine structure is a sum 
of the second-order perturbation corrections induced by the 
Breit-Pauli Hamiltonian and the expectation value of the ef- 
fective fine- structure Hamiltonian to this order, 
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Here, l/(E — H )' is the reduced Green function and 
is the spin-independent part of the Breit-Pauli Hamiltonian, 
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where we omitted a term with <5 3 (r) since it vanishes for the 
triplet states. It is noteworthy that in Eq. (5) we include the 
operator Hf s [and not just H^], thus accounting for the amm 
correction to the Breit-Pauli Hamiltonian. While this correc- 
tion is of order m a 7 , it is convenient to calculate it together 
with the m a 6 contribution because only simple changes in 
the prefactors are required. 

consists of 15 operators first derived by Douglas and 
Kroll (DK) [16] in the framework of the Salpeter equation. 



The finite nuclear mass correction correction to the m a 6 
contribution can be divided into the mass scaling, the mass 
polarization, and the operator parts. The mass scaling prefac- 
tor is (m r /M) 4 for the B 2 , B 3 , B 4 , and B 5 , (m r /M) 5 for the 
other Bi operators, (m r /M) 6 for the second-order corrections 
involving the first term in Eq. (6), and (m r /M) 5 for all other 
second-order corrections. The mass polarization effect is most 
easily accounted for by including the mass polarization op- 
erator into the zeroth-order Hamiltonian. The operator part 
comes from recoil corrections to \ and The 

recoil part of is given by Eq. (4). The spin-independent 
recoil part of the Breit-Pauli Hamiltonian is 
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Recoil corrections to the DK operators were studied by Zhang 
[15] and by K.P. and Sapirstein [17]. The result is given by 
the effective Hamiltonian 
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tonian and is the sum of the spin-dependent and the spin- 
independent parts, — Hf^+H^2. The effective Hamil- 
tonian to order m a 5 is given by the sum H ^ = + H^l , 
where the first part is the leading-order amm correction to the 
Breit-Pauli Hamiltonian defined by Eq. (3) and 
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The amm part of H ( 5 ) has already been included into the m a 6 
correction described in the previous section. The remaining 
part of the second-order perturbation correction will be de- 
noted as Eg and is given by 
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The effective operator H^p consists of two parts: (i) the 
exchange terms, in which photons are exchanged between the 
two electrons, and (ii) the radiative corrections, in which one 
or several photons are emitted and absorbed by the same elec- 
tron. They are calculated separately using different computa- 
tional methods in the following subsections. 



IV. DERIVATION OF THE ma 7 CONTRIBUTION 

In this section we present a detailed derivation of the m a 7 
contribution to the helium fine structure. The corresponding 
results have already been presented in Ref. [9]. The derivation 
is based on the dimensionally regularized NRQED [18]. The 
general idea is that, in the situation when all relevant electron 
momenta are much smaller than the electron mass, an approxi- 
mate QED Lagrangian can be used, obtained from the original 
full-QED Lagrangian by the Foldy-Wouthuysen (FW) trans- 
formation, as described in Appendix A. The standard FW 
transformation is generalized to the extended number of the 
space dimensions and also to account for the magnetic mo- 
ment anomaly of the electron a e . The regularization parame- 
ter e, related to the space dimension d = 3 — 2 e, plays the role 
of both an infrared and ultraviolet regulator and cancels out in 
the end of calculations. 

The fine structure contribution to order m a 7 (a 5 Ry) can 
be written as [12] 

where i?W denotes the effective Hamiltonian to order ma 1 , 
1/(Eq — Hq)' is the reduced Coulomb Green function, and 
El is the low energy contribution to be interpreted as the rel- 
ativistic correction to the Bethe logarithm. 

The second term in Eq. (33) that involves the second-order 
matrix element is the simplest. is the Breit-Pauli Hamil- 



A. Photon exchange part 

An important feature that leads to a considerable simplifica- 
tion of the calculation of the photon exchange part is the fact 
that the order being calculated is nonanalytic in a 2 . For exam- 

(5) 

pie, H^ is consists of the two terms only, which can be derived 
from the two-photon exchange scattering amplitude. A sim- 

(7) 

ilar statement holds for the photon exchange part of H^ s : if 

(7) 

Hf s is an effective Hamiltonian, it has to give the same scat- 
tering amplitude as in full QED. Due to the simple structure of 

(7) 

Hf s ', the effective interaction can be unambiguosly extracted 
from this amplitude. The scattering amplitude is usually much 
simpler to calculate, than corrections within effective field ap- 
proach, such as that used for the derivation of H ^ . An impor- 
tant point here is that only the two-photon exchange diagrams 

(7) 

contribute to H fs ', the absence of the three-body effects being 
a result of an internal cancellation. 
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FIG. 1 : The two-photon exchange scattering amplitude 
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So, we obtain the exchange contribution from the spin de- 
pendent part of the two-photon scattering amplitude, which 

is 
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where q = p[ — p\. There are three scales of the k integral 
that are responsible for the to a 7 corrections: to, to a, and 
to a 2 . Only the first two scales are accounted for in Eq. (36), 
whereas the third one, k <~ ma 2 , corresponds to a low-energy 
contribution and requires a separate treatment. Because of 
the dimensional regularization, the contribution of each en- 
ergy scale can be obtained separately. 

In our calculation, only the spin-dependent part of the scat- 
tering amplitude 8\H is needed. It is, however, not obvious 
what the spin dependent part is. In order to be consistent with 
the rest of the calculation, we employ the free FW transfor- 
mation S, 



i> + i 



^/2E p (E p + l) ' 
which for small momentum takes a simple form 

2 ■ 

This leads us to the following projection operators 

u(p')Qu(p) = TrQu(p) ® u(p') 

fTrQ(4i)(!^(£±i), 
l^TrQ(ai) (4^)4(^1), 



(37) 
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which identify the spin-independent and spin-dependent parts 
of the matrix element of the arbitrary operator Q, respectively, 
with ct 4 - 7 defined by 
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We now perform an expansion of the integrand of the scat- 
tering amplitude 8\H in Eq. (36) for two scales, k <~ m 
and k ~ ma. Assuming k ~ to and the external momenta 



p ^ ma and expanding the integrand in a, we obtain 
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A similar expansion with the assumption that k ^ ma leads 
to 
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The sum of <5i_ff[m] and 5\H[ma] is 
= a 2 |CTi(i,g)cr2(j,9) 
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where <r(j, g) = cr JJ q % , and g = y<p. 

The third scale k <~ too; 2 requires a more accurate treat- 
ment since any number of the electron-nucleus Coulomb pho- 
ton exchanges contribute to the same order. This low energy 
part can be represented in the Coulomb gauge as 



E L e = e 
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where the symbol S (. . .) stands for the first-order perturbation 
correction of the matrix element (. . .) by the (<i-dimension 
generalization of the) Breit-Pauli Hamiltonian H^\ which 
implies perturbations of the reference-state wave function <f>, 
the energy Eq, and the zeroth-order Hamiltonian Hq. 

The expression (45) involves the Coulomb Green function, 
which is not known for the arbitrary dimension. This problem 
is solved by splitting the integral over k into two parts, 



dk = 



dk 



dk , 



(46) 



with A = m (Z a) 2 A and A being a dimensionless cutoff pa- 
rameter. The two corresponding parts of Eq. (45) will be re- 
ferred to as and S 2 E. It is assumed that in these parts 
the expansion is performed first in the small e and next in the 
large A. The first part £ LE has a finite limit at d = 3. It 
will be evaluated in Sec. IV E together with other low-energy 
contributions. 

We now turn to the evaluation of the second part 62E. The 
spin-dependent part of the Breit-Pauli Hamiltonian in d di- 
mensions is the sum of the electron-electron part H ee and the 
electron-nucleus part H c n, 
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So, the contribution due to the photon exchange is given by 
the sum of the expectation values of 5\H and 82H and by the 
low-energy contribution £le- 

When calculating expectation values of effective operators 
between the triplet P states, some simplifications can be per- 
formed. The first one is that the expectation value of the Dirac 
5 function with both momenta on the right (left) hand side 
vanishes. The second one is that the expectation value of 0^ 
is equal to that of a 2 - As a result, the sum of 5\H and 8 2 H 
can be written as 
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B. Radiative corrections 
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Only the electron-electron part H ee contribute to S 2 E. Since 
k is much larger than H n — E , we expand the integrand to 
yield 
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The first term in the above expression is the second-order 
perturbation correction, which is already included into Eg, 
Eq. (35). The contribution of the second term is 
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After expanding this expression in e = (3 — d)/2 and then 
in a, the result is represented by the expectation value of the 



The m a 7 contribution induced by the radiative corrections 
is also split into the high and low energy parts. We argue that 
the high energy part can be accounted for by using the electro- 
magnetic form factors F\ and F 2 and the Uehling correction 
to the Coulomb potential Fy, 

Fv(-f-) = -^9 2 . (53) 

7T 15 

In principle, there are also corrections quadratic in the electro- 
magnetic field. However, one can demonstrate that such terms 
formed out of E, B, p, and a contribute only to higher orders 
and can be neglected. 

Corrections induced by the slope of the electromagnetic 
form factors are obtained by rederiving the Breit-Pauli Hamil- 
tonian, Eqs. (47) and (48), with employing modified electro- 
magnetic vertices. The resulting effective operator is 

S 3 H = nZa(F{+2F 2 ' + F{ / )t[a 1 (p'{,Pi) + ^(p' 2 \p2)} 
-7T a(2 F[ + 2 F' 2 + F{r )i \cf\ (p\ , Pl ) + a 2 (p' 2 , p 2 )} 
-2na(2F[ + F 2 ' + F[ r )i[ ( j l (p l 2 ,p 2 ) + ( j 2 (p[ , Pl )} 
+7T a(2 F[ + 2F 2 + Fy) u\ (j, q) a x (j, q) , (54) 

where p" is the momentum scattered off the Coulomb poten- 
tial of the nucleus. 
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The low-energy part of the radiative contribution is written 
in a form similar to Eq. (45), 



erators and requires a proper definition. We introduce the fol- 
lowing operators 



E L r = e z 
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Here, S denotes the first-order perturbation correction due 
to both parts of the Breit-Pauli Hamiltonian, the electron- 
electron part SH ee in Eq. (47) and the electron-nucleus part 
8H e M in Eq. (48). Introducing the splitting parameter A, we 
separate El a into two parts, £lr and S4E, which correspond 
to the first and the second term in Eq. (46), respectively. 

The evaluation of S4E is similar to that of the photon ex- 
change part. It yields S4E = (S4H), where the effective 
Hamiltonian is 



S 4 H 
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The total radiative correction is the sum of S 3 H, 84H, and 
the low-energy contribution £lr- The sum of 83H and 84H 
can be simplified further by using the symmetry 1 <-> 2, with 
the result 
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C. Q-operators 



It is convenient to consider the sum of Eqs. (52) and (57), 
Hq = He + Hn, as several logarithmic terms cancel out. 
The results is 
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For numerical calculations, we have to obtain the coordi- 
nate space representation of Hq, which involves singular op- 
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The coordinate-space representation of these operators is de- 
fined through their integrals with the arbitrary function / that 
is smooth at origin, 
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where we assume that r is expressed in atomic units. 

With these definitions, the coordinate-space representation 
of Hq is (in atomic units) 
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+ 3 In a I ipi x 5 (r) P i ■ u\ 
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-a — «Pi x — pi • (7i 

4 7T r^ 1 



(63) 



The above equation is written in atomic units because the 
definitions (61) and (62) are formulated in this unit system. 
Other formulas in the present paper are written in relativistic 
units. The A-dependent term in Eq. (63) cancels with the cor- 
responding contribution in Eq. (122). The logarithmic in a 
part of Hq agrees with the results of Refs. [10, 11]. 



7 



D. Anomalous magnetic moment correction to the ma G 
operators 

(7) 

The remaining contribution to ' is the amm correction 
to the spin-dependent m a 6 operators. It does not lead to any 
divergences and therefore can be calculated without any reg- 
ularization. The derivation of the amm part of is done 
with the help of the NRQED Hamiltonian obtained by the 
FW transformation of the Dirac Hamiltonian with the electron 
magnetic moment anomaly included. The resulting Hamil- 
tonian, with higher-order spin-independent terms omitted, is 
[19] 



Hfw 



7T 2 



eA°- e -{l + a e )a-B-^ 



+ 



- (1 + 2 a e ) [V • E + a ■ (E x tt - n x E)] 
8 

(3 + 4a e ) 2 



32 



{p , e E x p- a} 



(64) 



We use the opportunity to correct the misprint in Ref. [9] 
where the last term was typed with an incorrect prefactor. As 
demonstrated in Ref. [20], all spin-dependent operators to or- 
der m a 6 can be obtained from Hpw- The derivation of the 
amm correction to the m a 6 operators is very much similar. 

We start with the general expression for the one-photon ex- 
change amplitude between the electron a and the electron b, 



(k) 



+ 



f b (k)e 



i k ■ r a 



E - H - k° + ie 



-ik-rh 



i k-r b 



1 



E -H -k° + ie 



fa(-k)e 



(65) 



where G^ v is the photon propagator in the Coulomb gauge, 



G„„(fc) 



i 

-1 



fc 2 



v = 0, 
i,v = 3 



(66) 



<j> is an eigenstate of H and j% is the operator of the electro- 
magnetic current for particle a. In the following, we will con- 
sider separately the exchange by the Coulomb Goo and the 
transverse photons. The expression for the electromag- 
netic current is obtained from the Hamiltonian Hpw as a 
coefficient that multiplies the electromagnetic potential A^. 
The first terms of the nonrelativistic expansion of the current 
are 



f(k) = l + -^a-kxp--^P 
4 m 8m z 



for the j° component and 



p_ 

m 



- — a x k + . 
2m 



(67) 



(68) 



for the j component. 

The main part of the calculation is performed in the non- 
retardation approximation, which consists in setting k° = 
in the photon propagator G^ u (k) and in the current j(k); 
the retardation corrections are considered separately. Em- 
ploying the nonretardation approximation and the symmetry 
k° «-> — k°, the integration over k° is carried out as 



1 r dk° 

2 J 2ni 



1 11 

-AE- k° + ie + -AE + k" +ie\ ~ ~2 ' 

(69) 

The one-photon exchange amplitude in the nonretardation ap- 
proximation thus is 

f a {k)e i% ^-^f h {-k) <f>Y (70) 

To the leading order, the current does not depend on k and the 
k integration gives the coordinate-space representation of the 
photon propagator in the nonretardation approximation, 



GV(f) = 



Gjj, v (k) 




0, 



i, v 



(71) 



One easily recognizes that in the nonrelativistic limit Goo is 
the Coulomb interaction. This term is already included in Hq, 
which means that the nonrelativistic Coulomb interaction has 
to be excluded from the perturbative expansion. Next-order 
terms resulting from the expansion of j° and j lead to the 
Breit-Pauli Hamiltonian. 

We are interested in the expansion terms that yield effec- 
tive operators of order ma 6 a e . Their derivation is analogous 
to that of the m a 6 Hamiltonian in Ref. [20], the only differ- 
ence being that the corresponding amm prefactors should be 
retained. These prefactors will give us the required effective 
operator, which will be denoted as . It is expressed as a 
sum of various contributions 



(72) 



which are calculated in the following. 

SHi is the correction due to the last term in Hpw in Eq. 
(64). This term involves only A and its gradient, so the 
nonretardation approximation is valid here. SHi includes the 
Coulomb interaction between the electron and the nucleus and 
between the electrons. We denote by V the nonrelativistic in- 
teraction potential 



a a>b b 



(73) 



and by £ a the static electric field at the position of particle a 



e£ a = -\7 a V = -Z a'-^ + ^]a^ , 



(74) 



b=£a 



' ab 
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and write 5Hi as 

3 + 4a 



independent on a e . 

8H 2 is the correction to the Coulomb interaction between 
electrons which comes from the 5 th term in Hpw, namely 



SHi = 32m 4 e ° a ' (Pa e £a xp a + e£ a xp a plj 

a 

= a e y g" Pi ^3 x Vi ■ °i + - Pl x Pl • aij , 

(75) 

If the interaction of both electrons is modified by this term, 
where by = we denote that the equation is valid modulo terms the nonretardation approximation holds and Eq. (70) yields 



--(l + 2o e ) V-E + a-(Exp-pxE) 



(76) 



J 



811,. XX / ' /:;/ ' ^ (1 + 6 2 4 ae)2 ( ^ + 2 ; 7 ^ / 7 " x 1 ) < ' ' " ( " - - - ; ' ^ ' ^ 



a>fc 6 

fa e 2 (l + 2a e ) : 



16 



J d 3 k ^ 



Ak-r 



ai-PiX he 1 " r - cti • pi x fc -p- S 2 ■ k x p 2 



(77) 



3i a ^ /^^s^ _ 3 a ^ _ ^ x ^ _ a , _ _ . , _ _ . 

Y ^5 r X ( r '-P2)Pl ' <7i - - X (r X pi ■ <7 1 )p2 ■ 2 - ^3 (Pl ' CT 2 ) (f> 2 ' ^l) 



where by = we denote the equation that is valid modulo spin 
independent terms. To make the comparison with previous 
calculations more transparent, we transformed operators to the 
same form as in the original DK derivation. 

SH 3 is the relativistic correction to the transverse photon 
exchange. The first electron is coupled to A by the nonrela- 
tivistic term, 



— p-A- 
m 



2 m 



S- B. 



(78) 



and the second one, by the relativistic correction, i.e., the last 
but one term in Eq. (64) , 

I ({<? ■ B, it 2 } + a e {7? • B, 7? • S}) | [p 2 2p ■ A 

+ 2p- Ap 2 + S ■ Bp 2 +p 2 S ■ B 

+ a e (p-Bp-S + p-Bp-B)], (79) 



It is sufficient to calculate 6H3 in the nonretardation approx- 
imation, in which any correction can be simply obtained by 
replacing the magnetic field A by the static field A a , 



r ab r ab 



+ 



2r «'\ ' ab 

a(l + a e ) (B b x f ab ) % 



ab 



The result then is 



Pb 



(80) 



^3 = ^1 [ZplPa ■ Aa + 2pa ■ Aapl + pl&a -V a X Aa + &a -V a X Aa pl + CLefa ' (V a X Aa) Pa ' S a 



+Pa -B a p a - (V a X A a )) 



(81) 



dk e 
2 

= a e a 



[2p\pi -Ai+p\Vi xAi- Si +a e pi ■ Sipi ■ Vi x Ai] 



3 2 r ■ Si f ■ S 2 1 
4 Pl ' r> 



f _ 1 
r> ' P2 ~ 4 1 



1 2 r _ _ '6 2 r ■ ai r ■ a 2 l _ _ _ r _ 1 _ _ 1 3^ 

2 p i ^3 x Pi ' CT i + 4P1 ^5 2 Pl ' CTlPl X ^ ' P2 ~ l Pl ' 01 Pl ' a2 ^ + 4 Pl ■ CTi:Pl ■ r ^ ■ CT2 



- 1 3 - 
" ^ + 4 r 



(82) 



where by °= we denote the equation which is valid on the level of the expectation value of the operator on the triplet- 
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state wave functions. More explicitly, following Douglas and 
Kroll, we use the symmetry (1 <-> 2) of the wave function to 
replace terms involving <7 2 by terms with a\ . 

The effective operator SH4 originates from the coupling 



e 2 (l + 2a e ) _ -. 



a ■ (E x A - A x E) . 



(83) 



present in the fifth term in Eq. (64). The resulting correction 
is obtained by replacing the fields E and A by the static fields 
produced by the other electrons, with the result 



5Hi 



E 



e 2 (l + 2a e ) 



£ a X A a A a X &a 



DK,a c 



3 Za 2 _ _ _ 3« 2 , ^ , 

4 2rfr 3 4 r 6 

Za 2 _ _ _ Za 2 

^-3- n x p 2 • ax + TT -3- ? r x n • (7i r • p 2 

2 



2r 3 r 3 



Q 

"2r 



4 r x p 2 • <ii 



The effective operator <5i?5 comes from the coupling 



— A 2 . 



(84) 



(85) 



present in the first term of Eq. (64). Again, in the nonretarda- 
tion approximation, the field A a can be replaced by the static 
field produced by the other electrons, 

5H 5 = £ y ^ °= e ' A * ~ "°e ^ r x • ^ . (86) 



The effective operators SH 6 and <5i?7 represent the single- 
and the double-spin part of the retardation correction to the 
nonrelativistic single transverse photon exchange. To calcu- 
late them, we have to return to the general expression for the 
one-photon exchange amplitude, Eq. (65), and take the trans- 
verse part of the photon propagator, 



SE 



d 4 k 



1 



(2n) 4 i (fc0)2 -k 2 + t, 



§*3 



k l k 3 
k 2 



+ (a^b). 



E Q -H Q -k° + i e fh{ k)e ^ 



(87) 



We assume that the product f a (k) f b (—k) contains at most a 
single power of k°. This allows one to perform the k° integra- 
tion by encircling the only pole fc° = |fc| on the Rc(fc°) > 
complex half-plane and obtain 



SE 



d 3 k 



(2 7r) 3 2fc 
f a (k)e^ 
+(a <-> b) , 



k l k 3 



1 



Ep-Hp-k 



f b (-k)< 



(88) 



where k = \k\. The retardation expansion of the electron 
propagator yields 

1 1 H - E (H - E ) 2 

Ep-Hp-k k^ k 2 fc 3 ' y ' 

The first term here contributes to the Breit-Pauli Hamiltonian 
and the second term, to E^K Taking the current f in the 
nonrelativistic form, the third expansion term is 



k l k 3 



X [pb + 



1 + a e 
2 

1 + a e 



(? a x V a e^- r °(Ho-E Q ) 



(90) 



This is the most complicated term among the amm correc- 
tions, so we describe its evaluation in detail. The correction 
is split into the double spin part SE & and the single spin part 

SE 7 , 

SE = SE 6 + SE 7 . (91) 
The double spin part is 

v X fc) • & X ^ ( e <^. (Hp - Epf e- i% - F > 



(92) 



fc 4 

We use the commutation identity 

e itf a (jffo _ E ^2 e -ik-f b ^ + (fl ^ fe ) 

= ^[e^,[{Hp-Ep) 2 ,e-^ 

= -\([£M^*"]]) (93) 

to express this correction as the expectation value of the effec- 
tive operator SH 6 , 



k a(l + a e f 



32 



Pi 



2 i i r% r3 
P2 >°1°2 — 



Za 2 cti • rdi ■ f a 



Pi 



, , - - \ - 3 (cti • f ) (<7 2 • r) 

+ (r • (7l) C7 2 IT T \ 'P2 



(94) 



The single spin part is 



- EE 

a#6 b 



ie 2 \ . . /■ d 3 k 



x (e« k ^ (Hp - Ep) 2 e~* k ^ a a x k ■ p b 
-p a ■ db x he'**' (Hp - Ep) 2 e- itFb ) . (95) 
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With the help of the integral formula one obtains 

i 



= EE 



a (1 + a e ) 



a>b b 
-[pa,V] 



<J a X 



Tab PI 
Tab 2 



[V,p b ] 



Pi 



"a x > 
Tab 2 



■Pb 



2 - r ab 

Pbi°b X 

Tab 



~Pa 



Pi 



B Tab P$ 

r ab ' 2 



Z a _ f\ Z a r fi _ 



+"2 Pi -cti -p 2 xpi - -p? — (f -p 2 ) (rxpi) -cti ) . 



(97) 



The effective operator <5_ffg represents the retardation cor- 
rection to the single transverse photon exchange contribution, 
in which one vertex is nonrelativistic, Eq. (48), whereas the 
second comes from the fifth term in Eq. (64), 

- e{1+ 8 2ae) a-(Exp-pxE). (98) 
With the help of Eq. (88), one obtains 

a^b b 

x Te ( ( gi fe ? ° ^ a x ^ a + ^ a x °° e * * Pa ) 

X E -H -k (* - I * X e_i ^) + h ' C - ' 

(99) 



In the expansion of 1/ (Sq — -Ho ~ m Eq. (89) the first 
term vanishes because of the Hermitian conjugation and the 
second term contributes to order ma 6 a e . After commuting 
(H Q — E ) on the left, one obtains the effective operator 5H 8 , 



SH S = (l + 2a e ) E^^a • (4 x A a -A a x£ a ) + 1 -^ 



Aa-Pa'X S a +p a X a a - A a ,p\ 



DK,a c - TT 

= 5H 5 + a e 



3 i a 



--^-PiZP 2 x pi • cti - L -^-p\ ^ (r-p2) (fx pi) ■ B\ + 



ia o 1 



9a ,(MK^f) 
~^Pi Is 



, <?i ) (pi • ? 2 ) - - Pi ^ x p 2 • a l - - Pl — x pi • ai 



(100) 



Finally, after some rearrangement, the total amm correction where 
can be expressed as the expectation value of 17 operators, 



17 



Eh = {Hf}) =Y,m , 



(101) 



Hi 



Z 2 fl _ 



Ol . 



(102) 
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n 2 = 



3Z fi 



r 



cti (f -p 2 ) . 



(103) 



1 r 

Wl6 = -TPl • CTl pi X -p 2 , (H7) 

4 r d 



3Z f _ f\ 

rls = — ~ • CTl -g • CT2 , 

4 r 6 



(104) h 17 = ipi-CTi(~pi-CT 2 l + 3pi-f^-CT 2 ).(118) 



W 4 = -^-7 r x p 2 • cti , 

2 r 4 



<w 3 - - - - 

^ 5 = _ 4^ r '' Tir ' ,T2 ' 



1 f 

n& = T Pi -3 x Pi ■ VI > 

4 H 



(105) 



(106) 



(107) 



The operators above are intentionally written in a form very 
similar to that for the DK operators given by Eqs. (8)-(22). 
For most of Hi, there is an one-to-one correspondence with 
the DK operators, all the difference being the overall prefac- 
tors. There are only three exceptions. The first one is that our 
operator H$ cancels out in the DK calculation, while DK op- 
erator H$ cancels out in our calculations. The other two are 
related to the different spin structure of the next to last term in 
Eq. (64), which leads to the operators Hie and H17. 

E. Low-energy contribution 



1 r 

'Hi = t Pi r x P2 ' CTl 
4 r 6 



4r r\ 



H 8 = ~— 3 X P2 ' CTl , 



1 1 

H 9 = --p 1 -?r-p 2 rxp 1 -a 1 , (110) 

2 r 6 



3 i 

H w = ~ A f *r x (t • P2 ) Pi ' CTl , (111) 

4r b 



Hi 



3 _ 



8r 5 



H 



12 = • CT2P2 • CTl . 

8r d 



21 n 1 _ _ 

W13 = — Pi-rr-CTir--CT 2 , 
16 r b 



(108) 



(109) 



r x (r x p\ ■ 3i)p2 ■ ct 2 , (112) 



(113) 



(114) 



The low-energy part El comes from the photon momenta 
region k < A, A being the cutoff parameter introduced by 
Eq. (46). El is a sum of the low-energy contributions due 
to the photon exchange (£le in Section IV A), due to the ra- 
diative corrections (£lr in Section IV B), and an additional 
uv-finite contribution that was not considered in previous sec- 
tions. El can be conveniently derived from the low energy 
form of the electromagnetic interaction Hamiltonian [21], 



Hr 



E 



— e e -> 

-er a -E-- — a l a riB l ~- — ^a a -E a xp a 

4 m 



2m 



(119) 

This choice of the starting point for the derivation is more 
convenient than the general nonrelativistic Hamiltonian in Eq. 
(A9) since it makes transparent the high degree of cancellation 
between various terms. Specifically, the contributions of the 
second and the third terms in Hi cancel each other and only 
the first term contributes to El, 



El = z-a 



d 3 k 



o (2 7r) 3 2fc 



xS( <j) 



(n + r 2 ) 



E-H-k 



(n + r 2 ) 



,(120) 



H14 = p\ ~v ' CTl Pl * CT2 ■ 

8 r 6 



(115) 



where 5(. . .) denotes the correction to the matrix element 
(. . .) due to the Breit-Pauli Hamiltonian Using the rela- 



tion 



i , 1 



Wl5 = o P\ ~~ (r • S2P2 • CTl + (f • CTl) (p 2 ■ CT 2 ) 



Cti T • a 2 r-p 2 ) , 



i _ 1 

[ff,ri+r 2 ] = (pi+P2)+i — o [piXCTi+p 2 XCT 2 ,i?o-^o], 

J m Am 2 L 1 



(121) 



(116) i?L can be transformed to the following compact form, 
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Er = 



+- 



2a 
3^ 

iZ 2 a 3 



3tt 



(pi+p 2 )(H-E) In 
(a i 



2(H - E) 



(Zaf 

<? 2 ) 



(Pi 



In 



2(tf 



+ P2) 
S )" 



ln(2A)-Za 2 (ift x «5 3 (n ) ft • S x ) 



{Zaf 



n ^2 



(122) 



where the term with the second order matrix element has been 
dropped out, as all such terms are included in Eg. The term 
with ln(2 A) in Eq. (122) cancels out with the corresponding 
contribution in Hq in Eq. (63). 



V. NUMERICAL EVALUATION AND THE RESULTS 
A. Nonrelativistic wave function 

In order to obtain the nonrelativistic wave function, we use 
the technique by Korobov [22, 23], in which the spatial part 
of the triplet P states is represented as 

N 

i=i 

x [ri cxp(-ai n - ^ r 2 - 7» r) - (1 <-> 2)] , 

(123) 

where r = \r\ — f 2 \. Real nonlinear parameters a.i, pi, and 7$ 
are chosen quasirandomly from the intervals 

cti G [A U A 2 ], 
Pi € [Bi,B 2 ], 

7* G [Ci,C 2 ], (124) 

with the parameters Ai i2 , £1,2, and Ci^ being subjects of a 
variational optimization. In order to enforce the proper behav- 
ior of the wave function (123) in the limits n — > 00, r 2 — > 00, 
and r — > 00, the nonlinear parameters are subjected to the 
condition 

{ai+A,ai+7i,A+7i}> V^^. ( 125 ) 

where i?i is the ionization energy of the atom. In order to 
reproduce the behavior of the exact wave function for small 
values of n,2 and r, the variational parameters Ai i2 , #1,2. 
and Ci j2 are allowed to take negative values. To make the 
basis set more flexible, multiple sets of the variational param- 
eters Ai^, Bi t 2, and Ci,2 are introduced. Namely, the dou- 
ble basis set was used in this work for the determination of 
the nonrelativistic wave function, and the triple basis set was 
used in calculations of corrections to the Bethe logarithm and 
second-order corrections. 

The calculation of matrix elements of the nonrelativistic 
Hamiltonian is performed with the use of the simple formula 



for the master integral: 

is? W **'- 



- (a + 0)(p +-,)(-, + a y 1 

Integrals with any additional powers of 7^ in the numerator 
can be obtained by differentiating with respect to the corre- 
sponding parameter a, p, or 7. Matrix elements of relativistic 
corrections involve additional inverse powers of n, r 2 , and r. 
They can be obtained by integrating with respect to the corre- 
sponding parameter. In fact, all matrix elements required for 
the evaluation of the relativistic, QED, and the finite nuclear 
mass corrections can be expressed in terms of the rational, 
logarithmic, and dilogarithmic functions of a, (3, and 7. 

The procedure of generating the nonrelativistic wave func- 
tion looks now as follows. For the initial set of parameters 
A 12 , B 12 , and C\^, the nonlinear parameters ojj, Pi, and 7; 
with i = 1, . . . , N are distributed quasirandomly. Then, the 
N x N matrix of the nonrelativistic Hamiltonian Hq is eval- 
uated. The linear coefficients c; and the reference-state eigen- 
value Eq are determined by using the inverse iteration method, 
with the LDU decomposition employed for the inversion of 
the Hamiltonian matrix. Then the procedure is repeated for 
a different set of parameters Ai i2 , Bi^, and Ci i2 , looking 
for the minimum value of the energy Eq. The minimization 
problem is rather noisy, as the functional contains many local 
minima. So, the simplest simplex-like algorithms of mini- 
mization are probably the most appropriate for this task. Our 
calculations were performed in the quadruple, sixtuple, and 
octuple arithmetics, which were implemented in Fortran 95 
by libraries written by V. Korobov [24]. 

Our result for the nonrelativistic energy of helium for the 
triplet P state and the infinite nuclear mass is (in atomic units) 

E(2 3 P) = -2.133 164 190 779 283 205 146 96 ±? . (127) 

This value is the upper variational bound for the energy ob- 
tained with N = 6600 basis functions, and the uncertainty is 
the extrapolated lower bound. The value in Eq. (127) is by 
about 4 decimals more precise than the previously best result 
of Ref. [8]. 

B. Angular momentum algebra 

In the approach employed in the present investigation as 
well as in the previous studies by K.R and coauthors [12, 28], 
all the angular momentum algebra is performed in Cartesian 
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coordinates. Tensor product of the 3 P wave functions is rep- 
resented [37] in terms of the spatial wave functions and the 
spin operator s = (<?i + <r 2 )/2 as 

\ 3 P )( 3 P \ = K)0'l ^ - A , (128) 



1^| 3 P 2 ,to)< 3 P 2 ,to| = 



(129) 
(130) 



x ^2 s 2 5 ij - 3s* s J + 2 s J s 1 ' 

(131) 

where \ j) denotes the state with the Cartesian index j and the 
normalization of the spatial wave functions is 



(132) 



For the calculation of the second order matrix elements one 
needs formulae for the spin product 



s l s 3 s k = cP fe s l + - e jkl s l s l + - e m s 3 s l + - e ijl s k s l 



and for spin traces 



Trs* = 0, 

Tr s i s j = 2 6 ij , 

Tr s*Vs fe = ie ijk , 

Tr sVs fe s z = S^S kl + S jk S l1 



(133) 



(134) 
(135) 
(136) 
(137) 



Using these formulae, all matrix elements can be reduced to a 
form involving the spatial wave functions only. For example, 
matrix elements of the operators in the Breit-Pauli Hamilto- 
nian can be expressed as 

(Q-s)j = ie jk i(j\Q k \l)uj (138) 
(s-Q-^j = (j\Q 3l \l)vj (139) 

where Q is an arbitrary symmetric and traceless tensor (Q j7 = 
Q l1 and Q kk = 0) and 



u.j = (1,1/2,-1/2), 
vj = (-1,1/2,-1/10). 



(140) 
(141) 



for J = 0, 1, 2 respectively. 



C. Leading-order fine structure 

The dominant contribution to the fine structure comes from 
the spin-dependent part of the Breit-Pauli Hamiltonian H{ s 
given by Eq. (1). With including the nuclear recoil effect, the 



leading-order contribution to the fine structure in helium is 

E fs (J) = (H is )j 



/TO r 
V TO 



3 a 4 



-E! (1 + a e ) 2 vj + E 2 (l + 2 a e ) uj 



4 \ m 
+E 3 [ 1 + - a e j uj + — E 4 (1 + a e ) uj 



.(142) 



The corresponding results for the large and the small fine 
structure interval are 



v i 



V\1 



V TO / 



^(l + a e ) 2 + ^(l + 2a e ) 



Ez,, 4 m E 4 . 

+ T (l + - ae ) + --(l + a e ) 



(143) 



TO r \3 

a RooC 

TO / 



-^(l + «e) 2 + f (l + 2a e ) 



E 3,, 4 s m E 4 ' 
+ T (1+ 3 ae) + MT (1 + ae) 



where the constants Ei are given by (in atomic units) 



Ei — 
E3 — 
E 4 = 



2 Z eji 
-3ejfe» 
4 Z ejki ( j 



x V] 



x (Vi - V 2 



3 x (Vi+V 2 ; 
' 1 



(144) 

(145) 
(146) 
(147) 
i). (148) 



Our numerical results for the constants Ei are 

E x = 0.180 220 618 632 744(10), (149) 

E 2 = -0.277 401 358 712 829(10), (150) 

E 3 = 0.411 999 963 626 094(25), (151) 

E A = 0.241 945 125 695 21 (6). (152) 

These results are obtained with including the mass polariza- 
tion term into the zeroth-order Hamiltonian and thus contain 
effects of the second and higher orders in m/M. 



D. ma contribution 

The contribution of order to a 6 to the fine structure of he- 
lium is represented by Eq. (5). The most difficult part of its 
numerical evaluation is associated with the second-order con- 
tributions. First calculations of the second-order corrections 
to the helium fine structure were performed by Hambro [25] 
and Lewis and Serafino [26]. Two decades later, Yan and 
Drake [27] did these calculations to a much higher accuracy 
and demonstrated that the first results were much less accurate 
than it was claimed. The nuclear recoil and the amm effects 
were included into the second-order corrections in Ref. [8]. 
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An independent evaluation of the second-order corrections 
(including the amm part but not the recoil effect) was per- 
formed by K.P. and Sapirstein [28]. In the present work, we 
re-calculate all m a 6 corrections, with the intention to inde- 
pendently check the numerical convergence of the previous 
results and, more importantly, to check the nuclear recoil ef- 
fect on the second-order corrections, which was previously 
calculated only by Drake. 

The second-order corrections in Eq. (5) are of two kinds: 
the symmetric and the non-symmetric one. The symmet- 
ric contributions are the most numerous ones since they in- 
volve the 3 P, P, 3 D, 1 D, and 3 F intermediate states. The 
derivation of the calculational formulas for them is relatively 
straightforward along the lines presented in Sec. VB. 

The numerical evaluation of the symmetric second-order 
contributions was performed by employing the variational op- 
timization of the nonlinear parameters of the basis set for the 
Green function. Convergence of numerical results is rather 
slow for the 3 P intermediate states because of the singular 
character of the Breit interaction. The convergence can be im- 
proved by introducing singular functions into the basis set, as 
in Ref. [27]. We, however, prefer to exploit the flexibility of 
the basis set (123) and emulate the missing basis functions by 
using very large exponents. In order to effectively span large 
regions of nonlinear parameters, we used non-uniform distri- 
butions of the kind [29] 



regularized expression for the second-order correction, 



a t = Ax + (*: 



l)A 2 



(153) 



with a — 2 and 3, where the variable U has a uniform quasir- 
andom distribution over the interval (0, 1). 

The non-symmetric second-order contributions involve 
only the 3 P intermediate states. Care should be taken in the 
numerical calculation of this part due to the presence of singu- 
lar operators, namely the Dirac 5 and the p 4 operators. While 
a straightforward numerical evaluation is possible, a much 
better convergence is obtained by transforming the singular 
operators to a more regular form. 

The treatment of the second-order correction involving the 
8 operator is based on the global representation of the 5 func- 
tion introduced by Drachman [30], 



(0|47T(5 3 (ri)|n) 



-(/n + /o)+ V V„- -V„ 
att 2 



(154) 



4tt [<5 3 (n) + <5 3 (r 2 )] 



- + - ] 4/o 



1 



(Eq — Hq) 
1 



7 Hf s 



7 H b 



a=l,2 



ri r 2 ) ■"'(E -H y 

ri + r 2 ) Va (Eq-Hq)' 



- + -) H (s ) + (- + -) (H is ) . (156) 
n r 2 J I \n r 2 j 



A regularized expression for the second-order correction 
with the operator — 1/8 (p* + p|) can be derived by using 
the identity 



1 / — *2 i -*2\2 i 1 ->1 -1 
"o (Pi +^2) + 4P1P2 



(157) 



and employing the Schrodinger equation to transform the first 
term, as described in Ref. [8]. 

For the numerical evaluation of the non-symmetric contri- 
butions we used the set of nonlinear parameters obtained by 
merging two subsets, one obtained by the variational opti- 
mization of the symmetric second-order Breit correction and 
another, by an optimization of the symmetric correction for 
the model perturbation SV = 1/rf + 1 jr\. The model poten- 
tial SV corresponds to the most singular part of the regularized 
S operator. 

The numerical results for the m a 6 second-order correc- 
tions are presented in Table I. Our values for these corrections 
in the non-recoil limit (both with and without the amm part) 
agree well with the results by Drake [8]. In the recoil part 
of the second-order corrections, we observe some deviation 
from Drake's results. To localize the source of the discrep- 
ancy, we separate the recoil contributions into 3 parts, which 
are unduced by the mass scaling, the mass polarization, and 
the recoil operators. The corresponding contributions to the 
vol interval from the 3 P intermediate states are 3.26 — 3.66 — 
0.69 = —1.09 kHz, to be compared with Drake's values of 
3.25 — 3.66 — 0.06 = —0.47 kHz. Analogous contributions to 
the 1/12 interval are 1.17-1.87 + 0.91 = 0.21kHz, to be com- 
pared with Drake's values of 1.17-2.15+0.91 = -0.07kHz. 
The contributions to the i/q\ interval from the 1 P intermediate 
states are -4.54 - 7.68 + 2.93 = -9.29 kHz, to be compared 
with Drake's values of -4.52 - 7.66 + 1.87 = -10.31 kHz. 
The total difference between our results and those by Drake is 
rather small numerically and does not influence significantly 
the comparison of theory with the experimental data. 



where 



f k = E k - V + ^ Vi • V 2 , 



(155) 



E. Relativistic correction to the Bethe logarithm 



and V = -Z/n-Z/n+l/r. Noting that /„ = E n -E +.fo, 
one can cancel the E n — Eq factor with the denominator of the 
reduced Green function in the second-order matrix element. 
Using the completeness of the eigenfunctions, we obtain the 



The relativistic correction to the Bethe logarithm is given 
by Eq. (122). The two A-independent parts in the right-hand 
side of this equation will be referred to as the Eli and E^ 2 
corrections, respectively. We will start our discussion with 
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TABLE I: Second-order contributions to order m a for the large Ooi) and small (^12) fine-structure intervals in helium, including the recoil 
and the amm corrections. Units are kHz. 



State 




recoil 


tinim 




V12 
recoil 




3p 


-4894.29(2) 


-1.09 


-14.40 


-1569.62 


0.21 


-3.20 


1 P 


6595.64 


-9.29 


23.28 


-6595.64 


9.29 


-23.28 


3 D 


26.33 


-0.01 


0.07 


50.50 


-0.03 


0.16 


l D 








22.20 


0.05 


0.05 


3 F 








52.24 


-0.01 


0.24 


Sum 


1727.68(2) 


-10.39 


8.96 


-8040.32 


9.51 


-26.02 


Drake [8] 


1727.58(4) 


-10.81(4) 


8.95 


-8040.38(5) 


10.19(11) 


-26.02 



the simpler part E^2- For the numerical evaluation, it is con- 
venient to transform this correction to the equivalent form, 



E 



L2 



_a_ r 

3ti" Jo 



dk 



k 2 L 2 (k)-A 2 



(158) 



where 

L 2 (k) 



i((pi+ P2) x 



Ho + k-Ec 



(pi + P2) ■ s 



(159) 



and A 2 is the leading term of the large-fc asymptotic expansion 
of k 2 L2(k), which has the form 



k 2 L 2 (k)=A 2 - 



B, 



■ + ■ 



C lnfc D 



E 



k 3/2 f,2 f,2 f,5/2 ■ 

The two leading asymptotic constants are evaluated to be 



B 2 



-1 ^ (pi x <5 3 (n)pi • s+p 2 x S 3 (r 2 )p 2 ■ s) 



3V2 



to that for the Bethe logarithm [31]. In fact, each angular- 
momentum contribution to L 2 (k) has the same form as for 
the Bethe logarithm, the difference being only the prefactors 
coming from the angular-momentum algebra. (It is important 
that the difference in the prefactors leads to the disappearance 
of the fc _1 term in the large-fc asymptotics of L 2 (k).) In order 
to perform the integration over k in Eq. (163), one needs to 
know the function L 2 {k) for a wide region of k. As noted 
in Ref. [32], there is no need to perform the full variational 
optimization of the basis for each value of k. The idea is that, 
having got the optimized set of nonlinear parameters for the 
basis at k = k\ and k = k 2 , for all k in between one can 
use the basis obtained by merging together the two optimized 
sets. The asymptotic behavior of the integrand L 2 (k) for large 
k, together with its value at k = 0, £2(0) = — ((ri + f 2 ) x 
(pi +P2) -s), served as useful tests of the numerical procedure. 

The general evaluation scheme is as follows. First, we per- 
form a careful optimization of nonlinear basis-set parameters 
for several distinct scales of k: ki = 10*, with i = 1, . . . , i max 
and « max = 4. The optimization is carried out with incre- 
menting the size of the basis, until the prescribed accuracy is 
achieved. The size of the optimized basis employed in actual 
calculations varied from N = 600 for k\ = 10 to N — 1600 
(162) f° r t ' le wave an d ^4 — 10 4 ' yielding the numerical accu- 



(160) 



(161) 



Here we correct the overall sign in Eqs. (159) and (161) as 
compared to Ref. [12]. It is noteworthy that the k~ 2 asymp- 
totic behavior of L 2 (k) arises through an internal cancellation 
of the three angular-momentum contributions ( 3 S e , 3 P e , and 
3 D e , where e stands for the even parity), since each of them 
separately falls off as k^ 1 only. 

In order to accurately perform the integration in Eq. (158), 
we transform this expression to the following form 



E 



L2 



_-H.fr 

+ 

J K 



dk k 2 L 2 (k) 



k 2 L 2 (k)-A2 



B 2 

k 3/2 



2B 2 

A 2 k+—= 



(163) 



where k is a free parameter. 

In the numerical evaluation of E^ 2 , we exploit the fact that 
the integrand L 2 (k) obeys the variational principle, similarly 



racy of about 10 digits for L 2 (k). For each particular value of 
k < 10 4max , the calculational basis is obtained by merging the 
optimized bases for the two closest ki points, thus essentially 
doubling the number of the basis functions. According to our 
experience, such merging usually yields an additional digit of 
accuracy. 

The integral over k e [0, k] in Eq. (163) was calculated 
analytically, after performing the full diagonalization of the 
Hamiltonian matrix and using the spectral representation of 
the propagator. This allowed us to avoid problems associated 
with the pole on the real axis coming from the 2 3 S state. The 
parameter n was set to n = 10. The integral over k G [k, 00) 
was separated into two parts, k < 10 lmax and k > 10 Wx . The 
first part was evaluated by using the Gauss-Legendre quadra- 
tures, after the change of variables t = 1/k 2 . The second part 
was evaluated by fitting the integrand to the form 



,9 , B 2 lnfc /1\ 1 ( \ 

kL ^- A2 -k^ = ^{k) + ¥ w ^7=k 



(164) 
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where W(x) denotes a polynomial of x. For fitting, we used 
the function L 2 (k) stored on the interval k = 1, . . . , 100. The 
total number of fitting parameters in the above expression was 
about 9—11. The optimal form of the fitting function was 
selected by demanding it to reproduce the known asymptotic 
constants A 2 and B 2 for the function L 2 (k). The error due 
to the fitting procedure was estimated by comparing the in- 
tegration results for the fitted function and for the numerical 
integrand outside the fitting region, i.e., for k € [10 2 , 10 4 ]. 
Our results for the asymptotic constants A 2 and B 2 are 

A 2 = 0.120 944 339 354 433 (8) uj , (165) 
B 2 = -0.982 581 108(2) uj . (166) 

The final result for the Eli correction to the helium fine struc- 
ture is (in units m a 7 ) 



E L2 = 0.067 682 1(5) uj , 



(167) 



This is in reasonable agreement with the value obtained previ- 
ously in Refs. [12, 28], which is —0.06775(5) uj, except for 
the overall sign, which we correct here. 

We now turn to the evaluation of the Eli correction. In 
terms of the integral over the photon momentum, it is written 
as 



En = lim 

3 IT if— >oo 



,K 

\ dk k L\ (k) — Ai K — Bi InK 
Jo 



(168) 



where 



Li(k) 



2(H, 



(4) 



1 



fs (Eo-Ho) 

x (p\+pi)) + ((p{+pi) 



7(P\+Pl) 



1 



Ho + k-Eo 
1 



H + k-E 



<4 4) > - ^ 



1 



H + k-E 



(p\+pi) 



(169) 



and A\ and B\ are the leading terms of the large-fc asymptotic 
expansion of the integrand, 



, A B 1 C Dink E 

kL 1 (k) = A 1 + — + ^ m + + ¥ + 



with 

Ax 
Bi 



2 ( H, 



(4) 



1 



H, 



(4) 



(E — H 
1 

(E — H 



7 (ft +&)' 



1 ^Z[5 i (r l )+5 i {r 2 )] 



(170) 



(171) 



(Pi x $ 3 (ri)pi -s+p 2 x 8 3 (r 2 )p 2 ■ s) . 

(172) 

At k = 0, the integrand L\ can be evaluated analytically to 
yield Li(0) = 0. 



For the numerical evaluation, Eq. (168) is written in the 
form similar to Eq. (163), 



E, 



Ll 



+ 



dk 



fcLi(jfe)- Ai 



Si 




The main difference of the numerical evaluation of Eli from 
that of E L2 is that the integrand L 1 (k), contrary to L 2 (k), 
does not obey the variational principle [i.e., there is no func- 
tional whose minimum yields the exact value of L\(k)]. Be- 
cause of this, in evaluation of E L \ we have to use the varia- 
tional optimization results for the nonlinear basis-set parame- 
ters obtained for El2- This is a serious drawback since it is 
clear that the optimal set of parameters for the integrand L 2 (k) 
is not exactly optimal for L\(k), because of an additional sin- 
gularity introduced by the perturbing Hamiltonian H^f 1 . After 
some numerical experimenting, we found that this additional 
singularity can be well accounted for if the calculational basis 
for each k is not just doubled by merging two sets optimized 
for two scales k x and k 2 , but tripled, with the third part ob- 
tained from the second one by (alternatively) scaling the pa- 
rameters at and by a factor g = 10. This trick was inspired 
by the method described in Ref. [33]. 

With this modification, our numerical evaluation of E L \ 
was done similarly to that for Eli- Because of the tripling 
of the basis set, we used the optimized parameters with some- 
what smaller number of the basis functions but increased the 
high-energy cutoff parameter up to ks = 10 5 . Our results for 
the asymptotic constants A\ and B\ are 

A x = -0.028 038 047 8 (10) u,j + 0.054 037866 (4) Vj , 

(174) 

B x = -0.169 127 85 (20) uj + 0.146477 680 (2) Vj . 

(175) 

The final result for the Eli correction to the helium fine split- 
ting is (in units m a 7 ) 



E 



Ll 



-0.107 664 (6) uj + 0.118 4044 (4) vj . (176) 



The spin-orbit part of the above result is by about 25% 
larger than the previously reported value of Ref. [12] of 
— 0.0817(20)uj, whereas the spin-spin part is by about 10% 
larger that the previous value of 0.0959(4)vj. The reason for 
this deviation lies in the insufficient accuracy of the previous 
calculations. 



VI. SUMMARY AND DISCUSSION 

The summary of all contributions available for the fine 
structure of helium is given in Table II. Numerical results are 
presented for the large z/qi and the small v\ 2 intervals, defined 
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TABLE II: Summary of individual contributions to the helium fine 
structure. Units are kHz. 







" L 2. 


Z7"(4) 


29 618 415.54 


O OAT 71 *7 OO 

2 297 717.82 


E ( ' 


-1 556.97(2) 


-6 544.93 


£ (7) [log] 


82.59 


-10.09 


Subtotal 


29 616 944.16(2) 


2 291 162.80 


Drake [8] 


29 616 943.40(6) 


2 291 163.40(13) 


£<o[nolog] 


21.73 


7.42 


E H 


—4.21 


4.05 


£s[nolog] 


11.42 


— 1.21 


Eli 


—31.51 


—4.99 




4.01 


9.22 


S (7) [nolog] 


2.04 


14.48 


Higher orders 


±1.6 


±1.6 


Total theory 


29 616 946.20 ±1.6 


2 291 177.28 ± 1.6 


Experiment [7] 




2 291 177.53(35) 


Experiment [5] 


29 616 951.66(70) 


2 291 175.59(51) 


Experiment [4] 


29 616 950.9(9) 




Experiment [6] 


29 616 952.7(1.0) 


2 291 168.(11.) 



by 

14,1 - [P(2 3 P ) -P(2 3 P)]//i, (177) 
v 12 = [£7(2 3 Pi) - E(2 3 P 2 )] /h . (178) 

The parameters used in our calculations are: a -1 = 
137.035 999 679(94), cR^ = 3 289 841 960 361(22) kHz, 
and m/M = 1.370 933 555 70 x 10~ 4 . In the table, the 
correction E^ is given by Eqs. (143) and (144) and the cor- 
rection p( 6 ), by Eq. (5). P (7) [log] denotes the sum of the 
logarithmic parts of Es and Eq. The corrections Eq, Eh, 
and Es are given by Eqs. (63), (101), and (35), respectively. 
The complete listing of numerical results for individual terms 
contributing to Eq and Eh can be found in Ref. [9] and is 
not repeated here. The relativistic corrections to the Bethe 
logarithm E L1 and E L2 are given by Eqs. (168) and (158), 
respectively. 

The result for E^ in Table II is consistent with that of 
Ref. [9] after accounting for the newer value of the fine struc- 
ture constant. The result for E^ differs slightly from the 
corresponding value in Ref. [9], mainly because of the change 
in the recoil second-order Breit correction, which was previ- 
ously calculated only by Drake [8]. As can be seen from the 
table, different theoretical predictions that include contribu- 
tions up to order m a 7 log a and m 2 /M a 6 (entry "Subtotal") 
agree at a sub-kHz level with each other. 

The nonlogarithmic correction to order m a 7 has not been 
checked independently. The compilation of results presented 
for this correction in Ref. [8] is in part based on the derivation 
by Zhang [13-15], which was shown to be not entirely consis- 
tent [9], and in part includes calculational results by K.R and 
Sapirstein [12]. 

Table II shows that the calculational error of our results 
is almost negligible as compared to the experimental uncer- 
tainty. There is, however, a much larger theoretical error in- 
duced by the higher-order corrections. It was believed previ- 
ously [10, 12] that the higher-order ma 8 effects contribute 



well under the 1 kHz level. Particularly, the analysis pre- 
sented in Ref. [12] identified several m a 8 corrections that are 
enhanced by ln(Za) but nevertheless contribute only about 
0.1 kHz. In the present investigation, we found several non- 
logarithmic corrections that might contribute at the 1 kHz 
level. 

The first contribution comes from the mixing between the 
3 Pi and 1 Pi levels. The nonrelativistic 2 1 P i - 2 3 Pi energy 
difference of 61.3 x 10 6 MHz acquires the relativistic correc- 
tion of —17.1 x 10 3 MHz. If we consider the ma 6 second- 
order Breit correction to the energy of the 2 3 Pi state with the 
2 1 Pi intermediate states, the modification of the 2 : Pi — 2 3 Pi 
energy difference by relativistic effects alters the value of the 
correction by about 1 kHz. We thus estimate the theoretical 
uncertainty of the i/qi and z/12 fine structure intervals due to 
the mixing between the 3 Pi and 1 P\ levels as ±1.0 kHz. It 
should be mentioned, however, that all corrections due to the 
mixing cancel identically in the sum of the large and small 
intervals, v Q2 = ^01 + ^12- It is, therefore, likely that the theo- 
retical value for the interval v 02 is more accurate than that for 
the interval vq\ and v\ 2 separately. 

The largest identified m a 8 contribution to the interval z/02 
comes from the one-photon exchange diagram, which was 
evaluated to all orders in Za but to the leading order in 1/Z 
in Ref. [34]. The result obtained in that work for u a2 is 
0.1033 m a 8 Z 7 . Because of the Z 7 enhancement, the numer- 
ical contribution for helium is quite large, 13.1 kHz. This 
is, however, only the leading term of the \/Z expansion; the 
complete contribution for helium is going to be much smaller 
because of the screening. In order to estimate the screening 
effect, we compare the complete contribution to order m a 6 
for 1^02 in helium, which is —8.11 MHz, with the correspond- 
ing one-photon exchange term, which gives 86.4 MHz. The 
resulting estimate is ±1.2 kHz. 

The total theoretical error due to the higher-order effects 
specified in Table II for the z/ni and v\ 2 intervals is obtained by 
adding quadratically the two error estimates discussed above. 
We observe that for the small interval, the theoretical value 
agrees well with the experimental results, whereas for the 
large interval, a disagreement of about 3 standard deviations 
is present. It should be noted that the present theoretical un- 
certainties are much larger than those specified in previous 
investigations, the reason being that in most previous cases, 
the uncertainties represented the calculational errors only. 

Commenting on the situation when theory agrees with ex- 
periment for one fine-structure interval and disagrees for an- 
other, we have to state that we do not have any satisfactory 
explanation for it. All effects contributing to one interval con- 
tribute also to the other, both contributions being comparable 
in magnitude. We thus see no reason why a theoretical pre- 
diction for one interval should be significantly more accurate 
than for the other. Presuming that the experimental value for 
the z/ni interval is correct, we have to conclude that the ex- 
cellent agreement of our theoretical value for the v\ 2 interval 
with the latest measurement by Borbely et al. [7] is probably 
accidental. 

Finally, we present separately the theoretical result for the 
sum of the large and the small fine structure intervals in he- 
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hum, z/ 2(theo) = 31 908 123.5 (1.2) kHz. For the reason 
discussed above, its uncertainty is smaller than for the v 01 and 
v\2 intervals. The theoretical value disagrees with the exper- 
imental result of ^ 02 (exp) = 31 908 126.78 (94) kHz [5] by 
about 2 standard deviations. 

To conclude, we performed an evaluation of the helium fine 
structure that is complete to orders m a 7 and m 2 /M a 6 . Our 
results for the ma 4 , ma 5 , and m a 6 contributions agree with 
those reported in previous investigations at a sub-kHz level. 
The present evaluation of the relativistic corrections to the 
Bethe logarithm significantly improves upon the original cal- 
culation. The corresponding results reduce the previously re- 
ported discrepancy between the theoretical predictions and the 
experimental results. However, the remaining difference for 



the vqi interval is larger than the estimated contribution of the 
higher-order effects. This discrepancy needs to be resolved in 
order to make possible the determination of the fine structure 
constant by means of the helium spectroscopy. 
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APPENDIX A: FOLDY-WOUTHUYSEN 
TRANSFORMATION IN d-DIMENSIONS 

The Foldy-Wouthuysen (FW) transformation [35] is the 
nonrelativistic expansion of the Dirac Hamiltonian in an ex- 
ternal electromagnetic field. Following Ref. [36] we extend 
this transformation to the case where the dimension d of space 
is arbitrary. The Dirac Hamiltonian in an external electromag- 
netic field is 

H = a-Ti + l3m + eA a , (Al) 
where 7? = p — e A, 

<"*-(£ o)'"-{l-,)' ^ 

and 

{a i ,(T j } = 26 ij I. (A3) 
The FW transformation S [35] leads to a new Hamiltonian 

H FW = e lS (H-id t )e- tS , (A4) 

which decouples the upper and the lower component of the 
Dirac wave function up to a specified order in the 1 jm expan- 
sion. Here we calculate i?FW U P to terms contributing to the 
m a 6 correction to the energy. We use a convenient form of 
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the FW operator S, which can be written as 

1 



S 



13 (a ■ tt) 3 

I , . . „, P _ 



2 m im z 



J_e_ 

4 m 2 

i e 
3m 



i e 



24 m 3 
3 {(a.^) 2 ,a.S}^. 



[a • 7T, [a • 7r, a ■ E}] 



(A5) 



The FW Hamiltonian is expanded in a power series in S 

6 

H FW = Y, hU) + ■■■ 

3=0 



(A6) 



where 



= [iS.wW-ift], 

= j [i5 , ,W (j " 1) ] , for.? =2... 6, (A7) 

and higher order terms with j > 6 are neglected. The calcu- 
lation of nested commutators is rather tedious but the result is 
simply 



H FW = eA v + 



, (a-n) 2 (a-T?r , (<? • 7?) 



2m 



8 m 3 



[<7 • 7T, U ■ hi\ 



8 m 2 

i e 



16 m 5 



128 m 4 

% e 



16 m 3 

[ir • ir, [a • ir, [a • ir , a • E}]] 



16 m 4 



{(<J-n) 2 , [a-Tr,a-E]}. (A8) 



There is some arbitrariness in the operator S, which means 
that is not unique. The standard approach [35] relies on 
the subsequent use of several FW transformations and yields 
a result that agrees with the d = 3 limit of Eq. (A8) up to a 
transformation with an additional even operator. 

Our aim is to obtain a Hamiltonian suitable for calcula- 
tions of the m a 7 contributions to energy levels of an arbitrary 
light atom. In this case one can neglect the vector potential 
A in all terms having m 4 and m 5 in the denominator. Less 
obviously, one can also neglect terms with a ■ A3 ■ E and 
B 2 . This is because they are of second order in electromag- 
netic fields, which additionally contain derivatives, and thus 
contribute only to higher orders. After these simplifications, 
-Sfw takes the form 



FW 



2 m Am 8 m A 



V-E + a i: > {BP, 7T J }) 



8 m 2 



16 m 3 
3e 



16 m 3 



3 fa d > E ) 

e 

32m 4 l " ~ y ' y ' ' 128m 1 



'' {p 2 ,W 2 A°} + p6 



where 



64 m 4 



16 m 5 ' 



B 13 
E l 



&■ A* -d j A\ 

-vm° - 



b 2 ,b 2 ^ ]] 

(A9) 

(A10) 

(All) 
(A12) 



